Efficient Direct Tomography for Matrix Product States 
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In this note, we describe a method for reconstructing matrix product states from a small number of efficiently- 
implementable measurements. Our method is exponentially faster than standard tomography, and it can also be 
used to certify that the unknown state is an MPS. The basic idea is to use local unitary operations to measure in 
the Schmidt basis, giving direct access to the MPS representation. This compares favorably with recently and 
independently proposed methods that recover the MPS tensors by performing a variational minimization, which 
is computationally intractable in certain cases. Our method also has the advantage of recovering any MPS, while 
other approaches were limited to special classes of states that exclude important examples such as GHZ and W 
states. 



Matrix product states (MPS) [I-?] are a variational class of 
states that can be specified by a small number of parameters. 
Their importance in quantum many-body physics and quan- 
tum information science stems from the fact that they seem 
to capture the low energy physics of a wide range of one- 
dimensional systems. Indeed, White's density matrix renor- 
malization group [ ] can be seen as a variational method over 
this class of states [7-9]. 

Quantum state tomography [10] is a procedure that uses 
many copies of a physical system and performs a complete set 
of measurements on them in order to learn their quantum state. 
Because a generic quantum state of a system comprising n 
particles is described by a number of parameters growing ex- 
ponentially with n, tomography requires exponentially many 
copies of the system on which to perform exponentially many 
distinct measurements. In addition, the outcome of these mea- 
surements must be processed on a computer in order to in- 
fer the state of the system; this post-processing also requires 
an exponential amount of time. While this method is clearly 
not scalable, it plays a central role in benchmarking quantum 
information processors and also in establishing a direct and 
complete link between theoretical models and physical sys- 
tems. 

Matrix product states, on the other hand, are described by 
a number of parameters polynomial (typically linear) in n. 
Thus, if a physical system is believed to be in an MPS — for 
instance a low temperature spin chain or a quantum simulator 
meant to prepare such a state — it should in principle be pos- 
sible to learn its value from a small number of measurements. 
This would be a great advantage over standard quantum state 
tomography and would provide an essential tool for the exper- 
imental study of large scale MPS. 

In this note, we present a tomography protocol to learn an 
MPS using a small number (polynomial in ri) of efficiently im- 
plementable measurements. In addition, the post-processing 
required to infer the MPS from the measurement outcomes 
also scales polynomially in n. This contrasts with recent inde- 
pendent proposals for MPS tomography [11, 12] for which the 
post-processing needs to solve a problem which is in general 
intractable (NP-hard in some cases [13, 14]). Essentially, our 
method works by measuring the MPS tensors directly, while 



the other methods work by estimating the local density ma- 
trices of the state, then performing a constrained variational 
minimization over all MPS states. 

Our method has the additional advantage of recovering 
any MPS, while the other approaches were limited to spe- 
cial classes of states (injective MPS [15]) that exclude impor- 
tant examples such as GHZ and W states. As a drawback, 
our method requires quantum circuits of size polynomial in 
n, while the methods of [ 1 1 , 12] use quantum circuits whose 
size is logarithmic in the bond dimension \. (Note, however, 
that the modification proposed in [12] to handle GHZ and W 
states does require quantum circuits of size poly(n).) 

Direct MPS tomography In what follows, we describe our 
method for an open-boundary chain of n particles, where each 
particle is a e?-dimensional qudit. Generalization to the case of 
a periodic boundary ring is straightforward (with an increased 
measurement and post-processing complexity). Let be the 
matrix product state of this system, with bond dimension x- 
We are given many copies of \ip), and our goal is to determine 
its MPS representation. 

The key fact is that has Schmidt rank at most x> i- e > f° r 
any partition of the chain into two segments { 1 , . . . , I } and 
{I + 1, . . . , n}, the state can be written as 



x 



b i+i,... n) 



en 



where |<^'"'' ) are orthogonal states of the I first particles and 
|</>j +1 '"' n ) are orthogonal states on the last n — I particles. 
Now, look at the first k sites, choosing k such that d k ^ 1 > x- 
The reduced density matrix on these sites has rank at most x> 
so by performing a unitary operation, we can rotate the state 
into the subspace |0) ® (C d )® fc_1 . This disentangles the first 
qudit from the rest of the chain. Now we set aside the first 
qudit, look at sites 2, . . . , k + 1, and recurse. 

In this way, we get a sequence of unitary operations 
Ui, U~2, ■ ■ ■ U n -k+i, where U{ acts on sites i, . . . , i + k — 1. 
This sequence of operations transforms l^) into the state 
|0)®™~ |^?)> where \rf) is some pure state on the last k — 1 
sites. 
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Note that this is precisely the procedure for preparing an MPS 
state, running in reverse [16]. From this decomposition of \ip), 

H)=U 1 - 1 ---K\ +1 \0r n - k+1 \r 1 ), (2) 

one can write down the MPS tensors. We have 
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Jzi, . . . ,Z n ), (3) 



Zi,...,z„e{0,...,<2-1} 



where 



(z u ...,z n \Ur 1 ---U-l k+1 \0f n - k+1 \ V ). (4) 



As we will show below, this can be written as a product of 
operators acting on (C d )® fe_1 . 

In detail, our method is as follows. Let k = [log d x] + 1. 
For i = 1, . . . , n — k + 1, do the following steps: 

1. Prepare the state \ip') = f/,_i • • • Ui\i/j). 

2. Perform standard tomography on sites i through i + 
k — 1, to recover the reduced density matrix p' = 

Tvx,..„i-l;i +k ,...,nW)m. 

3. Let \(f>'j) (j = 1, . . . , rank(p')) be the eigenvectors of 
p' . Note that rank(p') < x < d k ~ 1 . Define additional 
vectors |<// } for j = rank(p') + 1, . . . , d k , such that the 
\4>'j) form an orthonormal basis for (C d )® fc . 

4. Define the unitary 

^EEi^'H^-^i- (5) 

a=0 j=l 

acting on sites i through i + k — 1. 

We can reconstruct the MPS representation of \i[>) as fol- 
lows: For i = 1, . . . , k — 1, and z € {0, . . . , d — 1}, define 
to be an operator on (C d )® fe_1 , which is given by 



I? = Ji.....<-i ® (|0><^|) ® J<+i fc - 



(6) 



T/ acts nontrivially on the i'th qudit, and acts as the identity 
on the other qudits. 

Fori = 1, . . . , n— k+ 1, and z 6 {0, . . . , d— 1}, we define 
V* to be an operator on (C d )® fe_1 , which is given by 

V? = ® (^^(JO)! ®I 2 ,..., fe ). (7) 



Here we view U~ x as an operator on (C d )® fe , and on the left 
we apply (z\ on the last qudit, and on the right we apply |0) 
on the first qudit. 

Then equation (4) can be rewritten as: 

c Zu ..., Zn = (or" 1 !? 1 ■ ■ -T£?V? ■ ■ ■ V:i k+1 \rj). (8) 

So T* and Vf are the desired MPS tensors, with boundary 
conditions given by (Ol®*^ 1 and |?/). Note that the bond di- 
mension is d , 1 , which is at most d ■ \- It is easy to check 
that this procedure uses poly(n, x) local unitary gates, and the 
classical postprocessing takes time poly(n, x)- 

Accuracy There are a few factors that limit the accuracy 
of MPS tomography in an actual experiment. In particular, the 
state of the system is generally not exactly an MPS, but merely 
close to one; and the measurements are noisy and restricted to 
finite precision. We will address each of these issues in turn. 

First, we slightly modify our tomography procedure. Note 
that the reduced density matrix on sites 1, . . . , I will actually 
be full-rank, though most of its probability mass will lie on a 
subspace of dimension at most x- So, each time we apply a 
unitary Ui, we also want to truncate the reduced state on sites 
i, . , . , i+k— 1 to a subspace of dimension d k ~ 1 . A simple way 
is to measure site i in the standard basis, and postselect on the 
0} outcome (which occurs with high probability). Note that 
we know explicitly how large is this truncation error, since we 
know the reduced density matrix on sites i, . . . ,i + k — 1. 

Clearly, if is close to a matrix product state, the trunca- 
tion errors will be small. More importantly, if \ip) is far from 
any MPS, the truncation errors will be large. Thus, we can 
do certified tomography, i.e., we can test whether an arbitrary 
state is close to an MPS. 

Now consider the effect of noisy measurements and finite 
precision. (Suppose, for simplicity, that the state is ex- 
actly an MPS.) Our estimate of the density matrix p' on sites 
i, . . . , i + k — 1 will be slightly wrong, and hence our esti- 
mate of the unitary Ui will be slightly wrong. In addition, the 
error in Ui can result in a truncation error when we measure 
site i and postselect on the |0) outcome. This truncation error 
affects the subsequent iterations of our procedure, where we 
estimate t/j+i, Ui+2 and so on. 

More concretely, let Ui be the unitary matrices for the true 
state \ip) and let Ui be the estimates produced by our proce- 
dure. Note that in order to represent the same MPS, these ma- 
trices need only to coincide up to some gauge transformation, 
which can be fixed arbitrarily. Equivalently, note that while it 
is crucial to obtain a good estimate of the d k ~ 1 -dimensional 
subspace on which p' is supported, it is less important to iden- 
tify the individual eigenvectors of p' exactly. Assume that we 
estimate this subspace within error e. Then, fixing the gauge, 
we have \\Ui — Ui\\ < i ■ e (where the linear dependence on 
i comes from the accumulation of truncation errors from pre- 
vious iterations). This implies that the reconstructed MPS \tp) 
will be close to the true state of the system, 



< n 2 e. 



(9) 
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Local measurements degeneracy The schemes for MPS 
tomography proposed in [11, 12] make use of local measure- 
ments on small intervals of the chain. The outcome of these 
measurements impose constraints on the global MPS. For in- 
stance, if one chooses these intervals to be of size 2k, the re- 
duced density matrix on any such interval has non-maximal 
rank due to the MPS structure. A minimization procedure is 
then employed to find the MPS that satisfied all these con- 
straints. 

We observe that this procedure will fail in cases where 
the MPS is not injective. This situation will occur for in- 
stance when the system is in the ground state of a degen- 
erate Hamiltonian. In those cases, the information gathered 
by any local measurement is not sufficient to recover the 
global state: the solution to the minimization problem is not 
unique, it is degenerate. A simple example is given by GHZ 
states of n qubits, i.e. \ip) = a|0)® n + e lv b\l)® n . Lo- 
cal measurements are insensitive to the relative phase ip, so 
the methods of [11, 12] fail to reconstruct this state. On the 
other hand, the final "boundary" state obtained at the nth 
step of our protocol will be (77) = a|0) + e lif \l), so the 
phase information is revealed during the final step of tomog- 
raphy. The same situation would occur with W states \tp) = 
^ (e ivi 1 100 . . . 0) + jf* |010 . . . 0) + . . . e^" |000 ... 1)). 

In [ 1 2], a fix was proposed to overcome this problem under 
the assumption that the degeneracy is independent of n. Their 
method is to identify the entire degenerate subspace and to 
perform further measurements inside this subspace to identify 
the state. We note that convolutional quantum error correcting 
codes [18, 19] are an important class of exponentially degen- 
erate MPS for which this approach would fail. Our method is 
suitable for this particular situation as well as any other MPS. 
In the case of convolutional codes, the unitaries Ui returned 
by our procedure would coincide with the encoding circuit of 
the code. 

Conclusion We have presented a very simple, direct 
method to perform quantum state tomography of matrix prod- 
uct states using a small number of measurements. In our 
scheme, there is no need to perform numerical optimization 
to recover the MPS from experimental data; this is its main 
advantage over recently proposed methods. It also has the ad- 



vantage of recovering any MPS, and is not limited to non- or 
constant-degenerate cases. Like other approaches, our scheme 
can be used to certify that the unknown state is indeed close 
to an MPS with given bond dimension. Finally, our scheme 
extends trivially to the case of mixed states for which the re- 
duced state on any interval has bounded rank. 
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